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The effect of confinement on the orientational structure of a nematic liquid crystal model has 
been investigated by using a version of density-functional theory. We have focused on the case of a 
nematic confined by opposing fiat surfaces, in slab geometry (slit pore), which favor planar molec- 
ular alignment (parallel to the surface) and homeotropic alignment (perpendicular to the surface), 
respectively. The spatial dependence of the tilt angle of the director with respect to the surface nor- 
mal has been studied, as well as the tensorial order parameter describing the molecular order around 
the director. For a pore of given width, we find that, for weak surface fields, the alignment of the 
nematic director is perpendicular to the surface in a region next to the surface favoring homeotropic 
alignment, and parallel along the rest of the pore, with a sharp interface separating these regions (S 
phase). For strong surface fields, the director is distorted uniformly, the tilt angle exhibiting a linear 
dependence with the distance normal to the surface (L phase) . Our calculations reveal the existence 
of a first-order transition between the two director configurations, which is driven by changes in 
the surface field strength, and also by changes in the pore width. In the latter case the transition 
occurs, for a given surface field, between the S phase for narrow pores and the L phase for wider 
pores. A link between the L-S transition and the anchoring transition observed for the semi-infinite 
case is proposed. 

PACS numbers: 61.30.Cz, 61.30. Pq, 61.30.Hn 



I. INTRODUCTION 



In the last years there has been a vast theoretical effort to understand the properties of liquid crystals in bulk 
phases. However, it is only recently that the interfacial phenomena presented by these complex fluids have attracted 
an increasing interest. The ability to manipulate the direction of the preferred molecular alignment, the director, 
by coupling to surfaces, is a property with a great practical interest. For this reason, there is a new effort trying 
to describe this rich surface phenomenology in the liquid crystals. These studies have mainly used phenomenological 
Landau-de Gennes approaches and aimed to describe, with more or less success, different phenomena such as subsurface 
deformations (ordering close to a surface), changes in the alingnment of the director due to the interaction with 
a solid substrate Q, nematic-isotropic transition ||, wetting ||], etc. On the other hand, microscopic-type theories 
claim to provide a more precise description of theses inhomogeneous problems since they take account of the structure 
at the molecular scale. In this framework, there is a recent literature describing attemps to answer such questions 
and it is supposed to provide more satisfactory theories in the interfacial phenomena. For example, by employing an 
Onsager-type theory R. van Roij et al. |7) investigated the problem of biaxiality and wetting for the Zwanzig model. 
P. I. C Teixeira has discussed, in a recent work the existence of subsurface deformations for the case of a confined 
nematic phase by using a generalized van der Waals theory. In previous papers, also employing this last approach, 
we studied surface phenomena in the case of a nematic in contact with a single wall Q as well as capillary effects in 
the case of a nematic in slab geometry with symmetric walls . In Ref. |l our attention was restricted to the study 
of a nematic in the presence of a solid interface. Model parameters where chosen such that the surface forced the 
molecules to lie perpendicular to the interface (homeotropic orientation) while in the nematic-isotropic interface the 
nematic director had a preference to lie parallel to the interface (planar orientation). In the model, surface parameters 
responsible for the order are the particle-wall interaction (ew) and the particle-particle interaction (ec) which played 
an important role in the interplay between anchoring and wetting transitions. A first-order anchoring transition 
between planar and homeotropic regimes was found; in a fi — T (chemical potential-temperature) representation (see 
Fig. 1) this transition appears as a line called anchoring line. This line approaches the nematic-isotropic transition 
line tangcntially (at constant fj,, eventually cuts this line at To), due to the existence of total wetting state at the 
coexistence. Later, the nematic phase was confined by symmetric walls, i.e both walls favored the same orientation 
(homeotropic) [ |Uj . Then, also the same anchoring line appeared and it was found to be very insensitive to the pore 
width except with respect to the localization of To that is shifted by capillary effects. These findings have a bearing 
on the new phenomenology which appears when the liquid crystal is confined. 

In the present paper, the theory is generalized to the case where the walls compete in molecular alignment. We 
believe that deeper insight into the surface-induced effects in liquid crystals can be obtained from a molecular approach 
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and, in this respect, density funtional theory appears to be a powerful tool. The confinement of the nematic phase is 
imposed by the presence of a boundary made up of two opposing flat surfaces, one favoring homeotropic orientation 
and the other favoring planar orientation. In these conditions, a spatial variation in the tilt angle between the director 
and the surface normal is found. The textures adopted by the nematic director depend very dramatically on the values 
of both surface field and pore width. For the case of a strong surface field ew and a wide pore this dependence is 
linear with respect to the coordinate normal to the interfaces except close to the surfaces. This configuration satisfies 
surface boundary conditions but entails an elastic energy which is minimum for a linear tilt configuration. For weaker 
surface fields and/or narrow pores a different texture occurs where the tilt angle adopts a steplike configuration. Our 
calculations also reveal the existence of biaxial behavior close to the surface favoring planar orientation as well as in a 
small region associated to the change of orientation in the steplike configuration. These biaxilities are direct surface 
effects due to the existence of different interfaces. 

In the next section we give a short account of the theoretical model used to describe the structure and thermody- 
namics of the confined liquid crystal. In section [II we show the different tilt configurations that result and argue 
that the phase transition between the linear-tilt and the steplike-tilt configurations can be related with the anchoring 
transition found in the semi-infinite problem. We also discuss how biaxiality is affected by the pore width. We end 
up with a section of conclusions. 



II. THE MODEL 



The theoretical model is a standard generalized van der Waals theory based on a perturbative expansion, using 
a hard-sphere (HS) fluid as reference system Ml. Details on the physical basis of the model and how to obtain its 
solutions numerically can be found elsewhere ]12||13| ] . Our starting point is the grand potential functional per unit 
system area A, Q[p]/A, whose functional minimum with respect to the one-particle distribution function, p{v, ft), 
which depends on both molecular positions r and orientations O, gives the equilibrium structure of the interface. 
This function, p(r, Cl) = p(z)f(z,fl), contains a mass distribution p(z) and an angular distribution /(z, These 
quantities vary locally with the distance from z = to z = H, H being the pore width. The expression for Q[p], in a 
mean field approximation is, 



tt[p]=F r [p] 



dvdv'dtldtl' p(r, Cl)p(r', ft) 



xv(r-r', &,&')- / / drdtlp(r,(l)[p-v w (r,d)}, 



(1) 



where p is the chemical potential and 

F r [p] = / drf hs (p(r)) + k B T I dvp{r) < ln(47r/(z, ft)) > 



is the reference system free energy. In the above expression fh s {p{^)) is the hard-sphere free energy density of a 
uniform fluid with a density equal to the local density at r and '< ... >' is an angular average. The attractive 
potential v contains anisotropic (dispersion) forces driving the liquid-crystalline behavior of the model material: 



v(r, n, n r ) = v A {r) + v B {r)p 2 {n • n') + v c {r) p 2 (n ■ f) + p 2 (n' • f 



(2) 



where f = r/r, and Va{t), vg{r), Vcif) are functions of the intermolecular center-of-mass distance r. Note that this 
potential is adequate to model uniaxial molecules with top-bottom symmetry. In this work we choose them to have a 
simple Yukawa form, i.e. Vi(r) = — e^exp (— A»(r — c))/r for r > a, and Uj(r) = otherwise, where a is the diameter 
of a hard sphere. 

The walls are modeled via the following potentials: 



v w (z,9) = -ewe- Xw{z " 7) P 2 (cos( 
which favors homeotropic anchoring, and 



(left wall, z = 0) 



(right wall, z = H) 



(3) 



(4) 



which favors planar alignment. The parameter ew is the surface strength of the walls and plays an important role in 
anchoring phenomena. 
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In order to describe the orientational structure of the fluid a tensorial order parameter is defined as, 

Qa ,= fdnf(z,Ci)( 3Cl ^- s ^ ) (5) 



where Sap is the Kronecker symbol and Q a is the a-component of CI. This quantity is a traceless symmetric tensor 
with eigenvalues Ai > A2 > A3. In this description, the isotropic state is defined by the condition Ai = A2 = A3 = 0. 
Uniaxial states occur when there is a twofold degenerated eigenvalue and then two cases can be distinguished: uniaxial 
nematic state (Ai > A2 = A3) and random planar (Ai = A2 > A3). In the first case, the molecules align preferently 
along the director, defined as the eigenvector corresponding to the largest eigenvalue (the uniaxial order parameter 
U = Ai). In the second case, the molecular orientations are uniformly distributed on a plane perpendicular to the 
eigenvector correponding to the lowest eigenvalue. 

For the case of non-degenerated eigenvalues, we can define the uniaxial (U) and biaxial (B) order parameters 
in terms of the eigenvalues of this tensor. The former is defined, as before, as the largest eigenvalue of the order 
parameter tensor (Ai) and the corresponding eigenvector is the local director of the fluid. The biaxial parameter is 
proportional to the difference of the two smallest eigenvalues. This quantity gives information of the amount of the 
orientational order on the plane perpendicular to the director. In this way, a secondary director can be defined as 
the eigenvector corresponding to the eigenvalue A2. By the symmetry of the problem, we have focused in the case 
where the director varies in the xz plane and is invariant with respect to a reflexion y — ► —y. In these conditions, 
the director is characterized by the tilt angle ip, defined as the angle formed between the director and the z axis. By 
convention, we define B = ±2(A 2 — A3)/3, with B > if the secondary director is on the xz plane and B < if it is 
along the y axis. 

For numerical reason, it is convenient to obtain the orientation distribution of the molecules referred to a laboratory 
reference system, described by three order parameters, 



T}(z) = j d<j> sin 6d6f(z,n)P 2 (cos 6) (6) 

v(z) = J d<f>sbi9d9f{z,Ci)sw20co8<f> (7) 

a(z) = J d<j) sin 6d6f(z,Cl) sin 2 6 cos 20 (8) 
These independent parameters can be related to the set {ip, U, B} by the following expressions, 

tan^(z) = V ^ Z) (9) 



(10) 



V (z) - <j{z)/2 + ^/[ V (z)-<j(z)/2f + V (zy 
U(z) = I (r,{z) + la(z) + 3 V / R*) - {l/2)a{z)] 2 + v{zf 

B{z) = \ (v(z) + \a{z) - v /[r,(z)-(l/2Wz)P + ^) 2 ) • (11) 

Numerical values for the potential parameters were taken as ca = 1 (which sets the temperature scale), es/^A — 
0.847 and ec/ € A = 0.75, The range parameters Ai are set, in units of a (throughout we choose this unit to set the 
length scale), to Ai = 2, 4, 1.75, i — A, B,C respectively, and A\y = 1- The model predicts a bulk phase diagram with 
vapor, isotropic liquid and nematic liquid coexisting at a triple point temperature Tniv flltl- 



III. RESULTS 



Confinement of simple liquids generally brings about capillary effects and associated capillary condensation phe- 
nomena whereby condensation occurs below the saturation point in bulk, as measured for example by the chemical 
potential. Similar phenomena must occur in nematic liquid crystals but with additional complicating factors due to 
the orientational order induced by the confining surfaces. 

In the following we use the chemical potential /1 as the external thermodynamic potential, in addition to the tem- 
perature T. Capillary condensation of the nematic from the isotropic phase occurs when, at fixed /i, the temperature 
Tni{H) at which the pore becomes filled with an oriented (nematic) phase is shifted with respect to the corresponding 
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temperature in bulk Tjvj(oo). We are thus seeking a non-zero value of AT = Tnj(H) — TV/(oo). Our calculations 
reveal that, contrary to the case of symmetric walls [[To) and, for large pore widths (H larger than 20a), no signifi- 
cant capillary effect exists, i.e. we find that the NI transition occurs at the same temperatures as in the bulk state. 
However, there are indications that for smaller pore widths this situation is not true longer (this aspect is being 
studied currently). Note though that condensation lines involving confined vapor states are expected to be affected 
by confinement; since these lines are not central to our argument we do not give any results in the following. 

Fig. 1 shows the bulk phase diagram, with vapor (V), isotropic (I) and nematic (N) phases and their corresponding 
phase transitions. We have superimposed the anchoring line, representing surface phase transitions between states 
with different director alignment for the semi-infinite problem (i.e. only one surface, inducing homeotropic alignment); 
these results have been obtained with a surface field ew = 0.53, for larger values this transition curve tends to approach 
the coexistence lines reducing then the region of planar states || . 

Now, we have calculated equilibrium states but with the sample confined between two opposing walls separated 
a distance H. There is a molecular ordering field due to the surfaces through the potential parameter ew- If the 
surface field strength is sufficiently high, one expects the director to adopt the configuration which satisfies the 
orientation favored by each surface, at least approximately. However, this necessarily demands that the director 
becomes distorted within the pore which entails a free energy cost due to the elastic contribution. In fact, the 
distribution of this inhomogeneity along the pore is found to depend sensitively on the values of pore width and 
surface field strength, and is expected to also depend on the thermodynamic parameters /x and T. 

Competition between these two energies results in the existence of two tilt configurations: L, where the tilt angle 
rotates uniformly along the pore and the resulting tilt-angle profile is linear in z; S, where the tilt angle adopts a 
step-like configuration, being constant throughout the sample except in a narrow region close to the z = surface 
where a crossover exists from 0° (homeotropic) to 90° (planar). 

Fig. 2 shows an L phase for H = 20er, ew = 0.7 and for thermodynamic conditions of temperature and chemical 
potential where, in the semi-infinite case, the nematic is clearly in the homeotropic state (/i = —3.7, T = 0.57). In 
the middle of the pore the tilt profile is almost a linear function of z, approaching smoothly the walls to their favored 
values. 

In the Fig. 3 the density, the uniaxial and biaxial order parameter profiles are plotted. The profiles present a smooth 
behavior, reaching a clear plateau in the middle of the pore with values close to the bulk values. The vanishing value 
of B in regions far from the interfaces is consistent with the uniaxial nature of the fluid particles. Close to walls the 
profile behaves differently depending on the favored orientation. Around the z = surface the U order parameter 
increases due to the strong surface interaction (despite of a decay in density). However, these deformations do not 
induce biaxiality (i.e. B = 0) due to the smooth character of the tilt profile along the pore. On the other hand, near 
the z = H surface, where the tilt angle is 90°, the U order parameter gets lower values and B < 0, that reflects the 
orientation of the molecules are mainly distributed on the xy plane, being the x axis the in-plane preferred direction. 
The decrease of the amount of order is due to surface potential that promotes a random planar distribution. 

If the surface field strength is reduced the director tilt configuration changes over to a step, i.e. a S phase. Here 
the tilt angle is uniform and equal to 90° except in a narrow homeotropic region close to the homeotropic wall (Fig. 
4). In general, this homeotropic region increases with the wall-particle potential range and with the strength of the 
surface field, so for the case of ew = one recover the configuration which corresponds to a nematic-hard walls 
interfaces (planar throughout the pore). Note that the S phase is spatially asymmetric, the tilt orientation being 
different from planar only in a small region close to the homeotropic wall. In fact, a corresponding S phase where 
the tilt were in a homeotropic configuration except in a very small region close to the planar wall is not observed 
in our calculations because of the effect due to the ordering field coming through the potential parameter ec- This 
field promotes molecular configurations where the director lies along a direction either perpendicular or parallel to 
the direction of inhomogeneity, depending on the sign of ec (ec > or ec < 0, respectively). In our case ec > 0, 
which means that the preferred director orientation is planar except in the region, close to z = 0, where the effect of 
the wall promoting homeotropic alignment is more pronounced. 

In the figure 5 we plot the density, uniaxial and biaxial order parameters profiles in the S phase obtained with 
ew — 0.30. Again the profiles only deviate from the bulk values close to the surfaces. The deviations around z = H 
has the same origin that in the L phase. However, close to z — this behavior has different nature. An analysis of 
the order parameteres profiles reveals that the system goes smoothly from an homeotropic to a planar configuration 
via a random planar state in the xz plane. This point (zi at fig. 5) corresponds to the jump of the tilt (origin of the 
step function profile). This mechanism is reponsible of the decay of the order and the biaxial behavior around this 
region. In fact, the existence of a crossover between homeotropic configuration and a planar configuration is closely 
linked to the capability of the fluid to have biaxial behavior in a small spatial range between both anchoring states. 
So, the roles that the z and x axes play as the local director and secondary director, respectively, for < z < z± are 
exchanged at z = z\. However, in order to make this mechanism feasible a dramatic depletion on the amount of order 
is needed. Moreover, the secondary director changes its alignment from the z axis to the y axis as z is increased (i.e. 
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B changes its sign), before that biaxiality virtually disappears in the middle of the pore. Biaxiality directly induced 
by the breaking of symmetry associated with the presence of the interfaces and, in the case of a free interface with 
planar director alignment, has been observed in simulations |l4| ]. 

In order to study in more detail the dependence of the nematic ordering on the pore width, we have calculated 
uniaxial order and biaxiality profiles for different values of H and for ew = 0.7 (Fig. 6). For pore widths H — 20, 30 
and 40, for which the L phase is the stable phase, the profiles exhibit inhomogeneity regions close to the surface 
favoring planar alignment, being almost constant in the rest of the pore. Note that biaxiality occurs only where the 
tilt angle is significantly close to 90°. An interesting feature is that the order parameter profiles converge around to 
2 = (resp. z — H) to the profiles obtained for the z = surface against a nematic in a homeotropic configuration 
(resp. the z = H surface against a nematic in a planar configuration). Furthermore, such semi-infinite profiles are the 
equilibrium ones for these values of ew-, A* and T. On the other hand, for H — 7.5, the stable phase is the S phase, 
and additional regions of inhomogeneity (and biaxiality) develop in the neighbourhood of the z = surface. 

The transition between the L and S phases involves a finite free-energy barrier, implying the existence of a thermo- 
dynamic first-order phase transition between the two phases. This is revealed by computing the relevant free energy 
for confined systems, i.e. the grand potential, as a function of the potential parameters and/or the thermodynamic 
parameters [i and T. In our case we have chosen to fix the latter and vary ew- Fig. 7 shows the behavior of the 
grand potential as a function of ew- Two branches, corresponding to the L and S phases, are shown. The point 
where the branches cross, e^, gives the transition point. Note the existence of metastable branches in both phases, 
which indicates the first-order nature of the phase transition. Tilt angle profiles along the pore corresponding to the 
coexisting L and S phases at ew — £\y = 0-554 for H = 20a are shown in Fig. 8. Note that the linear behavior in the 
L phase has an smaller range and its slope is also less than in the case presented in Fig. 2. On the other hand, the 
step in the S phase is located further from the z — surface than in the case presented in Fig. 4. 

As in any analysis of metastable states, some caution has to be exercised in analysing the occurrence of the different 
phases. In particular, the initial conditions required in our numerical minimisation scheme have to be chosen carefully. 
For example, the reality of the metastable L branch is shown by performing the following analysis: i) Equilibrium S 
states at low ew are obtained by conducting minimisation processes starting from linear and step-like profiles; these 
processes give the same equilibrium S states, ii) For higher, increasing values of ew, using step- like profiles or linear 
profiles as starting conditions give different final states beyond the transition point e^y, the final L states having lower 
grand potential energies than the final S states. This is the region of metastability of the S phase. Proceeding in the 
same way but from the L branch at high ew and decreasing ew provides the metastability region for the L phase. 

Fig. 9 shows the global phase diagram with respect to the parameters ew and H and for thermodynamic conditions 
T = 0.57 and // = —3.7. The line indicates the phase transition separating L and S phases. As the pore width is 
increased, the region of stability of the L phase is enlarged since the elastic energy associated with the deformation of 
the director along the pore (proportional to the integrated squared gradient of the tilt profile) decreases as opposed 
to the case of the S phase, for which this energy is essentially constant and seems to approach asymptotically to the 
value of ew — t-wi a ^ which the anchoring transition occurs for the semi-infinite case of the nematic against the z = 
wall. For narrow pores, we find a quick increase of e^r as the pore width is reduced, indicating that, in order for the 
L phase to become stable, a large surface field is required to compensate for the large elastic energy associated with 
a high tilt angle gradient. An interesting issue, which we have been unable to address due to technical reasons, is 
whether this first order phase transition may end up at a critical point as the pore width is reduced or even become 
second order through the appearance of a tricritical point. 

An analysis of the results quoted above shows that there is a relationship between the L-S transition and the 
anchoring transition observed in the semi-infinite problem. Let us consider the case of H 3> o. We have observed 
that the density and order parameter profiles deviate from the bulk values only in regions around of the surfaces, 
being the ranges of the inhomogeinities of the order of the molecular interactions. This feature is not surprising for 
the S phase, since in the middle of the pore the nematic is in a non-distorted tilt configuration (i.e. planar), so its 
tilt profile can be understood as two uncorrelated semi-infinite profiles linked via the bulk planar state. However, 
this fact is also true for the L phase even when the tilt distorsion is propagated through the pore (see discussion for 
the Fig. 6 above). This is the scenario that the phenomenological theory implicitly assumes: the excess free energy 
AJ7 = + pV (where p is the bulk pressure and V — AH is the volumen of the sample) can be split in two terms, 
an elastic term AJ7 e that comes from the distortions of the director field in a macroscopic scale, and surface terms, 
Af^"g and Af^™ 1 "^, that include the effect of the density and order parameter inhomogeneities in a microscopic 
scale close to the boundaries. The elastic term is given in terms of the elastic constants K\, K2, K3 and K24 by the 
Oseen-Frank free energy [|l5|,[l6| as: 

Atl e = y(V-n) 2 + y(n-Vxn) 2 + y(nxVxnf-(if 2 + K 24 )V-(nV • n + n x V x n) (12) 

where n(r) is the director field. Note that the splay-bend term is not included since K13 = in bulk jL7|. If we 
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restrict our study to the cases in which there is translational symmetry in the xy plane and the director is assumed 
to be in the xz plane (there is no twist), then n(r) = (cos(^(,z)), 0, sin(^(z)). In our model of uniaxial molecules with 
top-bottom symmetry K± = K 3 = —2-Kp\Ul J Q drr 4 VB(r), where pi, and C/j are the bulk density and uniaxial order 
parameter, respectively [ [l9|Jl7|| . So, the total excess free energy is given by: 

An = Anz r f mwm + Anz r MHW(H)) + k 1 j H dz ^, {z))2 (13) 
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where ip' = dip/dz. Note that the dependence of the surface terms on the contact values at the boundaries of the tilt 
angle and its derivative respect to z |lq| . Functional minimization of (|l3|) leads to an equilibrium tilt linear profile 
ip' = constant — 0(1/ H). The exact expression for the profile will depend on the surface terms of the free energy, 
so it is needed some information about them. The values of AQ S Z ^ can be obtained microscopically from the excess 
free energy corresponding to a case of the wall against a distorted nematic whose tilt profile goes asymptotically as 
ip(0) + ip'(0)z for a <C z <C H (this discussion is analogous for the z = H case). Since ^'(0) <C 1/cr the effect of 
the distortion far from the z = surface will be a perturbation respect to the case ip' = 0, that corresponds to the 
semi- infinite case studied in ||. In this Reference it was shown that if the wall promotes homeotropic anchoring, 
depending on the ratio between e\y and ec two local minima can appear for ip = 0° and ip = 90° with a free energy 
barrier between them. However, if the surface favors a planar anchoring only a global minimum at ip — 90° can 
exist. We can expand now the surface free energy in a series of both ip and ip' around (ip, ip') = (0°, 0) and (90°, 0). 
The special symmetry of our problem (rotational invariance around the z axis and the top-bottom symmetry of the 
molecules) implies that the linear term in ip' vanishes. From this analysis we can conclude that for H 3> o two 
different states can exist: a linear state that goes from a tilt close to 0° at z — to a value close to 90° at z = H , 
and a state with non-distorted planar tilt through the pore, and that correspond to the L and S states we have found, 
respectively. Furthermore, by using this phcnomcnologic approach a first-order transition can be predicted close to 
the anchoring line, inside the homeotropic region. The equation (p~3h, up to terms of the order of 0(1/H), is: 



Afi_ AnZ f (ip ,0) , A^ri(V2,0) ( K 



A A A 

where iJjq — or ir/2 for the L and S states, respectively. The transition occurs for 



8fAo:r/(V2,o)-A^r / (o,o)) 



which implies that AQ s z ™~q (ir/2, 0) > Afi*"Y (0, 0), i.e. the conditions should correspond to have an homeotropic 
configuration in the semi-infinite case, and the transition is first-order due to the free-energy barrier that exists 
between the two minima for the semi-infinite case. So, for H > H t an stable L configuration exists and for H < H t 
the S state is the stable one. From the Eq. ( |l5| ) it is clear that H t — > oo as the system approaches to the anchoring 
line (characterized by the condition AQ s z ™~q (tt/2, 0) = AO*"q (0, 0)), i.e. the L-S transition occurs for infinitely wide 
pores and, consequently, the anchoring line can be understood as the limit of the L-S transition for H — » oo. 

Strictly speaking, the scenario presented above is only true for H 3> a. However, our calculations reveal that the 
L-S transition observed by our Density Functional Theory (DFT) calculations is the continuation of the transition 
predicted for wide pores. First of all, the free energy of the S phase is expected to converge towards its limiting value 
given by equation (|l4| ) very fast, as it was observed in the symmetric walls case in fic|. For the L phase, we have 
compared the DFT results for ew = 0.7 with the limiting expression given by the Eq. (MI) (see Fig. 10). The values 
of /i and T are the same as in all our calculations, i.e. /i = —3.7 and T = 0.57. For such conditions, Ks = 0.615. 
On the other hand, calculations for the semi-infinite case show that the asymptotic value for H — > oo is given by 
AVL/AKbT = 1.841. It is clear that the values of Afl converge for relatively small values of H to the expression given 
by the phenomenologic theory, and in any case the behavior is qualitatively correct. ^From this analysis we conclude 
that the L-S transition is driven by the anchoring (surface) transition that occurs for the semi-infinite case of the 
nematic against the z = surface. 

IV. CONCLUSIONS 

In this paper we have employed a well-tested generalised van der Waals theory to investigate the effect of confinement 
on the thermodynamics and the microscopic structure of a nematic liquid crystal. The presence of two confining 
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surfaces favoring opposite orientations of the nematic director makes the order parameter, which is a tensorial quantity, 
exhibit all of its richness, giving rise to capillary effects not present in simple, non-orientable fluids. The competing 
conditions at the surfaces imply a deformed state of the director inside the pore. 

The main result of our study is that the system may adopt two possible states: one where the director changes 
orientation at an interface located close to the surface that favors homeotropic anchoring (depending on the wall- 
particle potential range), and one where the deformation takes place uniformly, giving rise to a basically linear tilt 
profile, with some degree of deformation close to the surfaces. The director may change from one configuration to 
another via a first-order phase transition. The interplay between these two phases may be understood in terms of the 
competition between the elastic energy associated with a deformed director state and the energy associated with the 
orienting surfaces. We have also shown that this transition is closely linked with the planar-homeotropic transition 
observed in the semi-infinite problem. In order to prove this aspect, a phenomenologic treatment has been employed. 
On the one hand, the theory shows how for big pores (H ^> a) the L-S first order transition can occur at threshold 
value of the pore width H t . On the other hand, the phcnomcnological approach also reveals that for H — > oo the L-S 
transition corresponds to the anchoring transition of the semi-infinite problem. We also conclude that the transition 
obtained by the DFT is the continuation to shorter scales of the transition predicted by the phenomenological theory, 
that in principle is only valid for much larger pore widths. 
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FIG. 1. Phase diagram containing bulk phase transitions and the anchoring transition for the case ew = 0.53 (they join at 
the temperature Tb). NI coexistence is denoted by the thick solid line. The dotted line is the anchoring line for the case of 
a nematic in the presence of a single surface which favors homeotropic orientation; this line divides the nematic phase into a 
region with homeotropic (_!_) and planar (||) nematic states. The thin solid line corresponds to the bulk transitions involving 
the vapor phase. 

FIG. 2. Tilt-angle profile for a confined nematic at T — 0.57 and p = —3.7 (L phase). Values for the surface field and pore 
width are ew = 0.7 and H = 20<r, respectively. 

FIG. 3. Density (p), uniaxial (U) and biaxial (B) order parameter profiles for the L phase, ew = 0.7 and values of T and p 
are given in text. 

FIG. 4. Tilt-angle profile for a confined nematic with surface field ew = 0.30 (S phase). The values of T, p and H are the 
same as in Fig. 2. 

FIG. 5. Density (p), uniaxial (U) and biaxial (B) order parameters profiles for the S phase, ew = 0.3 and the values of T 
and p are the same as in Fig. 2. 

FIG. 6. Uniaxial (U) and biaxial (B) order parameter profiles for different pore widths and ew = 0.7. The pore widths H 
are, from left to right: 7.5, 20, 30, and 40. The values of T and p are the same as in Fig. 2. 

FIG. 7. Thermodynamic grand potential as a function of the surface field ew for thermodynamic conditions and pore width 
as in Fig. 2. Branches corresponding to step-like (S) and linear (L) tilt-angle profiles are shown. The point at which these 
branches cross gives the transition point at e w — 0.554. The values of T, p, and H are the same as in Fig. 2. 



FIG. 8. Tilt-angle profiles for L and S phases at coexistence corresponding to pore width H — 20 and critical value 
e w = 0.554. The values of T and p are the same as in Fig. 2. 



FIG. 9. Phase diagram in the ew — H plane. The values of T and p are the same as in Fig. 2. 

FIG. 10. Comparison between the Density Functional Theory and the phenomenologic theory (expression in the text). The 
diamonds correspond to the excess free energy per area in units of KbT of a nematic confined by opposing walls in the L state 
with T = 0.57, p = —3.7 and ew = 0.7, calculated by using the Density Functional approach. The solid line is the surface 
tension given by the phenomenologic theory given by the Eq. (|14j). 
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